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We present contemporaneous, broadband, near-infrared spectroscopy (0.9 
- 2.45 jum) and H-band photometry of the black hole X-ray binary, XTE 
J1118+480. We determined the fractional dilution of the NIR ellipsoidal light 
curves of the donor star from other emission sources in the system by comparing 
the absorption features in the spectrum with field stars of known spectral type. 
We constrained the donor star spectral type to K7 V - Ml V and determined 
that the donor star contributed 54±27% of the H-band flux at the epoch of our 
observations. This result underscores the conclusion that the donor star cannot 
be assumed to be the only NIR emission source in quiescent X-ray binaries. The 
H-band light curve shows a double-humped asymmetric modulation with extra 
flux at orbital phase 0.75. The light curve was fit with a donor star model light 
curve, taking into account a constant second flux component based on the di- 
lution analysis. We also fit models that included emission from the donor star, 
a constant component from the accretion disk, and a phase-variable component 
from the bright spot where the mass accretion stream impacts the disk. These 
simple models with reasonable estimates for the component physical parameters 
can fully account for the observed light curve, including the extra emission at 
phase 0.75. From our fits, we constrained the binary inclination to 68° < % < 79°. 
This leads to a black hole mass of 6.9M < Mbh < 8.2M Q . Long-term varia- 
tions in the NIR light curve shape in XTE J1118+480 are similar to those seen 
in other X-ray binaries and demonstrate the presence of continued activity and 
variability in these systems even when in full quiescence. 

Subject headings: binaries : close - infrared : stars - stars : individual (XTE 
Jl 118+480) - stars : black hole 



Introduction 



XTE J1118+480 belongs to the class of transient low mass X-ray binaries (LMXBs) in 
which a late-type donor star transfers mass to its compact companion (either a black hole or 
a neutron star) through an accretion disk. It was discovered by the All-Sky Monitor (ASM) 
aboar d the Rossi X-ray Tim ing Explorer (RXTE) satellite when it went into outburst in 
2000 (jRemillard et al.ll2000i ). Because of its location along a sightline of low interstellar 
absorption (E[B-V] = 0.013), this system has been extensively s t udied both in outburst 



and in quiescence at multiple wavele ngths (e.g.. iHynes et al.ll2000l ; iMcClintock et al.l 12001 



Fender et al. 


2001; 


Chatv et al. 


2003) 



system had returned to quiescence established several dynamical properties of the binary, 
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including its orbital period (P or b = 0.17 days), donor star radial velocity semi-amplitude 
(Ka = 709 km s -1 ), and a mass function of the compact accretor of f(M) = 6 .27 ± .O4M 
dMcClintock et al.lboOlbl ; l\Vagner et al.lhoOli bonzalez Hernandez et al]l2008h . The mass 
function is the minimum mass of the compact object, establishing that XTE J1118+480 
contains a black hole. 

Accurate compact object masses are required to test models of format i on and evolution 

of black holes and neutron sta rs ( Brown et alll998 ; Fryer fc KalogeralboOl ; Nelemans fc van dan Heuval 



200 ll ; iBelczynski et al.l 120121 ). A black hole is characterized by its mass an d spin, and a 



reliab le estimate of the latter is intricately dependent on the former (e.g., iSteiner et al. 
20091 ) . These two parameters are critical in understanding space-time behavior near a black 
hole ( iMcClintock et al.l 1201 if ). X-ray binaries with black hole accretors h ave been used 
to test accretion disk and jet models (jYuan et al.l l2005t iMaitra et al.l 120091 ). to test mod- 
els of natal kicks and s pin-orbit misalignments between the black hole and the accretion 



disk (iFraeos et al.l 120091 ). and to predict observational tests of braneworld gravity models 



( Johannsenl |2009[ ). Observations of black holes in X-ray binaries are also important for 
understanding supermassive black holes in AGNs. Studies have revealed a "Fundamental 
Plane" of black hole accretion in which X-ray luminosity, radio luminosity, and black mass 
are related, both for Galactic black holes and their supermassive counterparts, implying that 
physical processes in AGNs can be illuminated by observat ions of similar processes in X-ray 



binar ies by an appropriate scaling of the black hole mass ( Merloni et al.l 120031 ; iFalcke et al. 
20041 ). 



Compact object mass determinations require values for the binary orbital period, the 
radial velocity semi-amplitude of the donor star, the mass ratio, and the inclination. The 
latter is usually determined by modeling the ellipsoidal modulation of the light curve caused 
by the Roche lobe-filling donor star. The light curve observations are often performed in 
the near- infrared (NIR), where the late- type donor star dominates. Many investigators 
have assumed a negligible contribution from non-stellar sources (e.g., the accretion disk 
and/or a jet outflow) at NIR wavelengths. However this assumption has been questioned 
based on fits to the NIR spec tra with donor star templates and on variability in the NIR 



light curves during quiescence (IFroning et al.ll2007l ; iReynolds et al.ll2008l ; ICantrell et al.ll2008 



20101 ). Unaccounted-for extra flux will lower the derived binary inclination from its true 
value and consequently overestimate the compact object mass, so accurate determinations 
of the relative contributions of the donor star and other emitters in the NIR are necessary 
to obtain correct mass values. In a recent study of the effects of system atic errors on the 
derived masses of black holes in X-ray transients, iKreidberg et al.l ( 120121 ) showed that the 
assumption of zero non-stellar light biases derived black hole masses to higher than their true 
values, often by significant amounts. They also cautioned that the non-stellar component 
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often varies over time and with orbital phase, requiring careful modeling of the observed 
light curve and contemporaneous acquisition of light curve and spectroscopic data to obtain 
unbiased compact object masses. 



Gelino et al.l (120061 ) obtained the most comprehensive light curve data set for XTE 
J1118+480, covering the B, V, R, J, H and K wavebands. They simultaneously modeled all 
the light curves to derive an inclination of 68lJ'o degrees and a black hole mass of Mbh — 
8.53±0.6 M . However, they assumed the non-stellar contribution in the NIR was negligible, 
<8% dilution of the donor star light. Their inclination value was lower than the values 
found by other investigators (ranging from 71° to 82°) using a variety of analysis methods, 



2001 




Zurita et al. 


2002b: 



Changes in the shape of the NIR li ght c urve over short time periods — the J-band light 



curve t aken by iMikolaiewska et al.l (120051 ) shows uneven peaks, unlike the symmetric ones 
seen by I Gelino et al.l ( 120061 ) five months earlier — suggests that variable, non-stellar emission 



is in fact present in XTE J1118+480 in the NIR. 

Motivated by the goal of determining an accurate black hole mass in XTE J1118+480, 
we obtained contemporaneous light curve and spectral data of this system at NIR wave- 
lengths to establish the non-stellar dilution from spectroscopy at the time of acquisition of 
light curve so that true binary inclination of the system can be obtained. In the following 
sections, we estimate the veiling caused by the non-stellar components at NIR wavelengths 
by examining the broadband spectral energy distribution and by measuring the equivalent 
widths of absorption lines in the photosphere of the donor star compared to field stars of 
known spectral type. We account for the extra NIR light when fitting ellipsoidal models 
to the observed light curve data in order to obtain a robust value of binary inclination and 
consequently determine the black hole mass. 



2. Observations 

We obtained contemporaneous spectroscopic and photometric observations of XTE 
J1118+480 using the the Gemini Near-Infrared Spectrograph (GNIRS) at Gemini-North and 
the Near-Infrared Camera & Fabry-Perot Spectrometer (NICFPS) on the 3.5-m telescope at 
Apache Point Observatory (APO). Photometry and spectroscopy were initially scheduled to 
be acquired on 2011 April 2 and 3 at both locations. We observed successfully at both sites 
on April 2. Additional light curve data were obtained on April 3 but due to bad weather at 
Mauna Kea, we did not complete the spectral observations until April 12. Table [1] lists the 
observations, exposure times, and orbital phase coverage for the program. 
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2.1. Spectroscopy 



We observed XTE Jl 118+480 using GNIRS at Gemini-North felias et alJ20o3 l. GNIRS 
was configured in the cross-dispersed mode with the 31.7 line mm -1 grating and the 0.30" 
wide slit, yielding R ~1700 while effectively covering the near-infrared range from 0.9 - 2.5 
/im. We obtained a total on-source observing time of 4.2 hrs. Data were taken in ABBA 
pairs at two positions along the 7.0"-long slit with a 310 sec exposure time at each position. 
The slit was oriented along the mean parallactic angle throughout the observations. An AO 
V type telluric star was observed hourly with the same configuration. From two nights of 
data we were able to cover all binary phases at least once, except phases <fi = 0.10 - 0.35 and 
= 0.72 - 0.76, which were covered twice. The raw data were reduced using the GNIRS 
tools within the Gemini-IRAF package, version vl.ll. The raw images were first treated for 
fixed pattern noise using the cleanirlll tool. The other data reduction steps consisted of flat- 
fielding, sky subtraction and wavelength calibration. Owing to the faintness of the target, 
we could not confidently extract the spectra from the individual sky-subtracted exposures. 
Instead, we combined every three exposures (adjacent in phase and equivalent to 5% of one 
full orbit) to perform reliable spectral extraction. The spectra extracted in this manner 
were minimally affected by orbital smearing effects when compared with the instrumental 
resolution of 177 km s" 1 . Flux calibration and telluric correction were performed wit h the 
xtellcor package within Spextool (version 3.4) (jCushing et al.ll2004t IVacca et al.ll2003l ). We 
shifted each of the extract e d spe ctra to the rest frame of the donor star using the orbital 
ephemeris of ICalvelo et al.l ( 120091 ). median-combined each order separately, and merged all 
the orders to generate the time-averaged spectrum of XTE J1118+480, which is shown 
in Figure [TJ The spectrum has been boxcar-smoothed by 2 pixels, corresponding to one 
resolution element. The spect rum in Figured] was corrected for interstellar absorption using 
Ay=0.066 ( iGelino et al.ll2006l ). We did not attempt to quantify the flux calibration accuracy 
since the spectral analysis of XTE Jll 18+480 does not require absolute flux values. However, 
we find that the mean H-band flux of the spectrum agrees with the photometry within ~ 



/o. 



2.2. Photometry 



We obtained H-band light curve data with NICFPS at APO flHeartv et al.l 120051 ). We 
used individual exposure times of 20 seconds for the target. A standard star for flux calibra- 
tion was observed every 1.5 hrs. We had excellent seeing conditions (~ 0.6") on the night of 



1 http: / /staff .gemini.edu / astephens / niri/clcanir / 
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April 2, but on April 3 the seeing conditions went from moderate (~ 0.8") to poor (> 1.0") 
after the first two hours of observing and the data were discarded. The binary orbital phases 
between <fi = 0.1-0.2 and <fi = 0.33-0.5 were not covered at all, but all other phases were 
covered more than once. To account for the dark current in the NICFPS exposures, dark 
frames matching the exposure times of the target and standard star were obtained at the 
beginning of each observing run and subtracted from the target /standard star exposures. 
Finally, the images were flat-fielded with sky-flats. We obtained five point dither images of 
the field, each offset by 20", and constructed a sky-flat by median-combining these exposures. 
To obtain sufficient signal to noise in our reduced data images, we median combined every 
three exposures (equivalent to ~ 60 seconds) to construct final reduced images. Aperture 
photometry was performed on the combined images using the IRAF phot package. For flux 
calibration, we used the stars in the field of XTE J1118+480 selected from the Two Mi- 
cron Sky Survey (2MASS) catalog ue as well as the ARNICA (Arcetri NICMOS3 camera) 
near-infrared standard star AS- 11 (IHunt et al.lll998l ). 



3. Analysis 

3.1. The Donor Star Spectral Type and Fractional Contribution to the NIR 

Spectrum 

The spectrum of XTE Jl 118+480 displays narrow absorption lines of neutral metals, 
including transitions of Al I, Na I, Mg I, Fe I, and Si I that are believed to originate in the 
photosphere of the donor star. We also detect broad emission lines of H I and He I from 
the accretion disk. The J-, H- and K-band spectra are shown in Figures [21 [3] and HI The 
error bars shown were calculated for each resolution element (2 pixels) by fitting straight 
lines through several continuum-dominated regions and obtaining the rms scatter about the 
fit. The continuum-dominated regions were selected after an inspection of similar regions in 
K5 V - Ml V template stars. We obtained signal-to- noise (S/N) estimates of ~ 18 in the 
H-band and ~ 15 in the J- and K-bands, respectively. 

In the J-band we detect a few absorption features in the 1.175 - 1.320 jum range, mainly 
blends of Mg I, Si I, Fe I, and Al I, and broad emission features of H I and He I. In the H-band, 
we observe narrow atomic features of Mg I, Al I, and Si I. In the K-band, the most prominent 
features used for spectral classification are the 12 CO bands, Na I, and Ca I, with Na I. The 
position of these lines are marked in Figure HI There is also a weak emission feature of H I 
at 2.16 ixm. Most of the absorption features are detected at low confidence or undetected 
in K. We do not detect any CO features in the K-band. Enhanced N V and depleted C IV 
and O V emission lines have been seen in the UV spectrum of XTE J1118+480, suggesting 
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that the accreting material has been CNO processed (IHaswell et al.ll2002l ). Therefore, our 
non-detection of the CO-bands in the K-band is unsurprising. In the H-band, we apparently 
detect the 12 CO (6,3) feature at 1.619 /mi, but given the relative oscillator strengths of 
the K- and H-band lines (the second overtone, Au = 3, bands of CO near the 1.1619 /im 
have oscillator strengths that are approximately 100 times smaller than the first overtone 
bands near 2.29 /im) and the presence of several other species contributing features in this 
wavelength region (including Ca, Fe, Ni, Si, and OH), we consider the apparent H-band CO 
feature spurious. 

The spectral type of the donor star in XTE J1118+480 has be e n broadly classified 



as K5 V - Ml V (IGonzalez Hernandez et all 120081: iGelino et all 120061 : iWagner et all 12001 
Frontera et aJboOlllMcClintock et al.ll200lh . To obtain a precise black hole mass, we want to 



constrain the spectral type of the donor star further and to estimate the relative contributions 
of NIR light from the donor star and other sources in the system. Towards that end, we 
examined the broadband spectral energy distribution (SED) of XTE J1118+480 compared 
to field stars. Figure [5] shows the spectrum of XTE J1118+480 compared to the spectrum 
of a K5 V star (in red) and an Ml V star (in green). The broad H-band bump (between 
1.5—1.7 /im) that is seen in the spectra of K-stars increases in amplitude wi th decreasing 



effect ive temperature and is attributed to the H opacity minimum at 1.6 /mi ( IRayner et al. 



20091 ). In XTE J1118+480, the size of this feature is consistent with a donor star spectral 
type later than K5 V. 

Earlier studies assumed that the donor star is the sole source of NIR emission in XTE 
J1118+480. If we follow that assumption and scale the field star spectra to match the 
observed flux in K (at 2.23 /im), we find that the template field star spectrum unphysically 
exceeds the observed spectrum at shorter wavelengths (>7% in both J and H) for K7 V or 
earlier spectral types. Thus, some dilution of the donor star flux is required if the spectral 
type is K7 V or earlier. Based solely on the broad spectral shape, an MI V or later donor 
spectral type could be responsible for most or all of the NIR flux, at least for wavelengths 

> 1.5/iHl. 

To constrain the donor spectral type further and estimate the fraction of NIR flux 
originating from the donor star (donor fractio n, /) in XTE J1118+480 , we followed the 
procedure outlined in iKhargharia et al.l ( 120101 ) and iFroning et al.l ( 120071 ). Specifically, we 
fit normalized template spectra of known spectral type to the normalized spectrum of XTE 
J1118+480. Dilution of the X-ray binary spectrum by non-stellar emission sources will have 
the effect of decreasing the equivalent widths of the lines in the observed spectrum compared 
to those in the templates. By scaling the templates to fit the observed line strengths, we can 
determine the donor star fractional contribution to the NIR spectrum of XTE J1118+480. 
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To normalize the spectra, we fit a spline function to the continuum, which was then 
divided out. The continuum points were selected by eye. The normalized spectrum of 
the template star was scaled by a fraction / that was varied from to 1 in steps of 0.01. 
The normalized, scaled template was subtracted from the normalized spectrum of XTE 
J1118+480 and the residual was computed. The value of / that minimized the deviation 
between the residual and the mean of the residual determined the best fit. For the template 
spectra, we used stars of spectral type K5 V (HD36003), K7 V (HD237903), M0/M0.5 V 
(HD209290) and Ml V (HD42581) (There is a disagreement in literature about the spectral 



type of HD209290 and hence we refer to it as M0/M0.5 V star lRavner et alJl2009 



2010l ). The template star spectra were obtained from the IRTF spectral librarj 



Koen et al. 



The fit regions that were investigated in the J-, H- and K-bands for the spectral types 
between K5 V - M1V are listed in Table [2] along with the best fit donor fractions, /, and 
uncertainties on those values. In the J-band, our fits were restricted to the region between 
1.20 - 1.314 (Urn within which we detected features containing blends of Mg I, Si I, Fe I and 
Al I. In the H-band, we fit features of Mg I, Al I and Si I over the wavelength range 1.475 
- 1.73 /xm. We only show fits in the K-band to the Na I feature, although the detection of 
this line is marginal. 

Since the S/N in our spectrum is low and the use of y 2 st atistics for the spect ral fitting 



are affected by systematic uncertainties (see the discussion in lFroning et al.ll2007l ). we have 
adopted a new technique to properly evaluate the noise associated with the spectrum of 
XTE Jl 118+480 and its contribution to the donor fractions. The procedure for this involves 
"scrambling" the XTE J1118+480 spectrum and fitting the template star spectra to many 
such randomly "scrambled" spectra in order to place robust error estimates on the best fit 
donor fractions. A detailed explanation can be found in the Appendix. The uncertainties 
calculated in this manner are shown in Table [2J 

We also examined if the calculated donor fractions were consistent with the overall shape 
of the broadband spectral energy distribution of XTE Jl 118+480. To answer that question, 
we scaled the flux of each field star by the average donor fraction in the H-band (near 1.6 
fim) and computed the value by which the average flux of the scaled template spectrum 
drops below that of XTE J1118+480 in both the J- and K-bands. Less emphasis was given 
on obtaining a K-band match due to the fact that we only had the Na I feature to compare 
with and this feature is at about the same level as the noise in the K-band. Depending on 
whether the corresponding drop in flux in the J-band was consistent with the average J-band 
donor fraction (within the bounds allowed by propagating the errors to calculate the average 



' http : //irtf web . if a.hawaii . edu/~spex/IRTF_Spectral_Library 
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donor fraction) computed from Table we could further constrain the donor spectral type. 
When a K5 V star was scaled to 43% of the XTE J1118+480 flux at the center of the H-band, 
it resulted in an average flux drop of 41% in the J-band (between 1.1 - 1.3 /xm) and 38% in 
the K-band (near 2.20 /jm), inconsistent with the average donor fractions expected from the 
fits to lines in those bands. A K7 V matches the drop in J-band flux but not in the K-band, 
while MO. 5 V/Ml V matches both J and K-bands when their H-band fluxes were scaled to 
the corresponding average H-band donor fraction obtained from Table [2j This suggests that 
K7 - Ml V is the most likely range of donor spectral types in XTE J1118+480. This is also 
consistent with the shape of the broad H-band bump near 1.6 /im discussed above. 

Thus, we conclude that the spectral type best describing the donor star in XTE J1118+480 
lies between K7 V - Ml V. We estimate the donor fraction by averaging the best fits to mul- 
tiple lines using the K7 - Ml V templates and propagating the corresponding uncertainties 
from Table [21 This leads to a H-band donor contribution of / = 0.50 ± 0.32 in the spectrum 
of XTE J1118+480 at the epoch of our observations. The upper panel of Figure [6] shows the 
normalized spectrum of XTE J1118+480 plotted (in black) over the normalized spectrum of 
a K7 V star (in red) while the lower panel shows the same comparison but with the K7 V 
star scaled to match the best donor fraction of 50%. 



3.2. Modeling the Light Curve to Obtain the Inclination 



Our spectroscopic observations of XTE J1118+480 were obtained in conjunction with 
contemporaneous light curve data. We phase-binned t he H-band observations using the 
orbital ephemeris from I Gonzalez Hernandez et al.l ( 120081 ). Figure |S] shows the H-band light 
curve after the data was combined into orbital phase bins of size A0 = 0.03. Each datum 
in the figure represents the mean of the points in the specific phase bin. The error bars are 
the rms scatter of the points about the mean in each bin. In cases where we did not have 
enough data points (<3) to place a reliable error estimate on the magnitude, the error bar 
for that datum was changed to the value of the largest error bar for the other light curve 
points. The orbital phases represent the standard convention wherein phase 0.0 is inferior 
conjunction of the donor star. 

The H-band light curve in Figure |8] shows a departure from the conventional ellipsoidal 
modulation expected from a Roche- lobe filling donor star: the maximum at phase = 0.75 
is higher than the maximum at = 0.25. Similar asymmetric m odulations were detected i n 
the quiesc ent J-band light cur ve of XTE J1118+480 obtained by lMikolajewska et al.l ( 120051 ). 
H owe ver . iGelino et al. obta ined NIR light curves of XTE J1118+480 four months 

prior to iMikolajewska et al.l ( 120051 ) and found no asymmetries in their NIR light curves. 
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Typically, this observed asymmetry in the light curve maxima is attr ibuted to the emission 



from a bright spot in the accretion d isk (IFroning fc Robinson! l200ll ) or from a dark spot 



on the donor star (jGelino et al.ll200ll ). Evidence for the presence of a hot spot in XTE 



during 



( 2004 ) found no evidence of a hotspot in their Doppler tomograms. ICalvelo et al. ( 20091 ) 



pot, but 


Torres et al. 


ICalvelo et al. 


( 


20od) 



speculated that variations in the mass transfer rate from the donor star during quiescence 
may cause the bright spot to be intermittent. 

In the following sections, we present two simple models for fitting the observed light 
curve of XTE J1118+480: a) a model incorporating a donor star with constant extra flux 
from nonstellar sources (e.g., the accretion disk and/or a jet); and b) a model incorporating 
a donor star and an accretion disk with a bright spot. We modeled the H-band light curve 
of XTE J1118 
in 



Zh ang et al 



480 u sing an updated version of the light curve synthesis code first presented 
(119861) that has been used to model the light curves in several other LMXB's 



(IFroning fc Robinson! l2001t iKhargharia et al.ll2010l ). The code accounts for the geometry of 
the binary system and then calculates the light curve by computing the temperature and 
intensity distribution across the surface of the Roche-lobe filling donor star and an accretion 
disk with a bright spot. It also takes into account gravity and limb darkening for each 
component of the system. The best fit light curve is obtained by minimizing the chi-squared 
value between the synthetic and the observed light curves. Below we discuss each model 
that was fit to the observed light curve of XTE Jl 118+480. 



3.2.1. 



Modeling the H-band light curve with a donor star and constant non- stellar flux 



We ran models to evaluate the effect of adding constant (i.e., present over the full binary 
orbit) non-stellar flux on the inclination of the binary. For modeling purposes, the donor star 
parameters were set as follows: T e ff = 4000 K, to reflect the average temperature from the 
spectral type range determined abo ve; a gravity darkening coefficient of 0.08, assuming that 
the donor has a connecti ve envelope (lLucylll967l ; ISarnalll989l ) ; and limb- darkening coefficient s 
from Claret et al. ( 1995 ). We adopted the m ass ratio value obtained by Calvelo et al. ( 2009h. 



Torres et al. 



which agre es with past e s timates obtained by lGonzalez Hernandez et al.l ( 120081 ) 
(120041 ) and IZurita et al.l (j2002bl ). To avoid the extra flux at <fi = 0.75, we modeled the light 
curve between phases <fi =0.0 - 0.50 only. By allowing the constant non-stellar fraction 
(1 — /) to vary within the limits obtained from spectroscopy, we modeled the light curve 
data to derive the corresponding best fit binary inclinations. Using this method, we find the 



binary inclination to lie between 68° - 89°. Figure [9] shows the best fit light curve 



68°) 
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obtained by modeling the observed light curve with a donor star and for a disk li ght fraction of 
fdixk. = 0.18 (xt = 0.51). We varied the mass ratio within the limits specified in lCalvelo et al. 



(120091 ) and the results were unchanged. From this, we conclude that accounting for the non- 
stellar flux within its uncertainty bounds obtained from spectroscopy places a lower limit on 
the binary inclination of i > 68°. 



3.2.2. Modeling the H-band light curve with a donor star and an accretion disk with a 

bright spot 

The H-band light curve of XTE J1118+480 shows asymmetry in the maxima of the 
peaks at phases 0.25 and 0.75, which is often attributed to emission from a bright spot on 
the accretion disk. Accordingly, we also fit the light curve with a model that includes a donor 
star and a cool opaque accretion disk with a bright spot on its rim. By this method, we do 
not attempt to constrain the physical properties associated with either the disk or the bright 
spot. (Indeed, there is no a priori reason to believe that the non- variable emission comes 
from the accretion disk rather than partially or wholly from a persistent jet.) However, we 
can derive a range of possible inclinations by modeling this system with reasonable estimates 
for the parameters associated with the accretion disk and bright spot. 

The donor star parameters were set at the same values used in the earlier models. Vary- 
ing the donor star parameters had negligible effect (< 1° change) on the derived inclination 
of the system and therefore we keep them fixed for these set of models and only vary the 
parameters associated with the accretion disk and the bright spot. For the accretion disk, 
we adopted an inner radius of 0.001-R^ (where represents the distance to the inner La- 
grangian poi nt) and an the oute r dis k radius of 0. 75 Rr„ ; these numbers are taken from the 



work done by I Wren et al.l (120011 ) and lCalvelo et al.l (120091 ). We set the flare half-angle of the 



disk to a small value of 1°. The bright spot was added to the rim of the disk and spans an an- 
gle in the azimuthal direction whose position <p spo t and width A<fr spot are variable parameters 
in the light curve fitting model. The other input parameters are the temperature of the disk 
{Tdisk) and the bright spot (T spot ). Both the accretion disk and the bright spot were assumed 
to emit as black bodies with single temperatures and linear limb-darkening coefficients. We 
varied T disk between 2500 - 4500 K, T spot from 5000 - 20000 K, <f) spot from 40° - 120°, and 



A&spn t from 5° - 15°. The choice of T d i S k was motivated by the work done in lReynolds et al. 



( 120081 ). who found that thermal emission from a cool outer accretion disk (T^ s fe ~ 2000 - 
4000 K) could be used to model the excess NIR emission in the multi-wavelength SED of 
several XRBs (including XTE J1118+480). 

Using these parameters, we fit models to the full H-band light curve. In Figure dOl 
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we show the best fit light curve (xt= 0.57), with an inclination of 78°, Tdisk = 3000 K, 
T sp ot = 12000 K, <j) S p t = 85° and A<p spot = 5°. From an investigation of our fits over the 
entire range of input parameters, we find that all fits have xt values between 0.60 - 0.80 
and produced lower/higher T disk values for correspondingly lower/higher T spot values; e.g., 
a model with Tdisk = 2500 K and T spot = 10000 K for an inclination of % = 74° gave a 
comparable fit (xl =0.60) to a model with Tdisk = 4000 K and T spot = 13000 K for an 
inclination for % = 81° (xt = 0.63). 

We do not attempt to set constraints on either the disk or the bright spot parameters 
based on the light curve models with the lowest xt value, except to note that a simple 
model incorporating a donor star and an accretion disk with a bright spot on the rim can 
account for the extra flux at <fi = 0.75. More important, we can set an upper limit on 
the binary inclination through the absence of eclipse features in the light curve. Extensive 
multiwavelengt h observations of XTE J1118+480 in ou tburst found no evidence of eclipses 



in this system ( jUemura et al.l l2000t IWood et al.l 120001 ) . From our modeling, we find that 



eclipse features start to emerge in the model light curves at i > 80° for an accretion disk of 
outer radius 0.75Rl- Hence, to maintain consistency with the non-detection of eclipses in 
XTE J1118+480, we find an upper limit to the binary inclination of < 79°. 

In modeling asymmetric light curve peaks in X-ray binaries, some authors have at- 
tributed the as ymmetry to a dark spot on the donor star rather than a bright spot on the 



accretion disk (jGelino et al.ll200ll ). In that case, our assumption in the previous section 
that modeling orbital phases <fi =0.0 - 0.50 is preferable to modeling the full light curve is 
exactly backwards and we would have underestimated the true amplitude of the ellipsoidal 
modulation (as we would then be modeling the hump artificially depressed by the stellar 
spot). If so, our lower limit to the inclination of > 68° is conservative but still encompasses 
the higher inclination that would result from fitting the larger hump in the light curve. 

In summary, by modeling the observed light curve with a constant non-stellar frac- 
tion obtained from spectroscopy, we found a lower limit to the binary inclination in XTE 
J1118+480 as > 68°. Additionally, an upper limit on the inclination is dictated by the ab- 
sence of eclipses, as < 80°. Hence, we have constrained the binary inclination to lie between 
68° < i < 79°. 



4. Discussion 



In Table [3], we summarize the current best estimates for the physical parameters in- 
cluding the black hole mass for XTE J1118+480. By performing spectroscopy of XTE 
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J 1 1 18+480, we find that the donor star contributes a fraction / = 0.50 ± 0.32 of the ri- 
band flux. By varying the constant non-stellar fraction within the limits obtained from the 
spectroscopy while simultaneously imposing the absence of eclipse condition, we obtained 
a binary inclination of 68° < i < 79 °. Usi ng this inclination range, combined with the 



orbital parameters from ICalvelo et al.l (120091 ). we find the mass of the black hole in XTE 
J1118+480 to be 6.9 < { Mrh /M^) < 8.2. Pr e vious analyses have found inclinations ranging 
from % = 60° to i = 8 2° jFrontera et al.l 200ll: IZurita et all 12002 bi llMikolaiewska et al.l 12005 



Khruzina et al.ll2005l ; iGelino et al.ll2006l ). Unfortunately, because of the relatively low S/N 
of our spectrum, we were only able to constrain the donor star fraction broadly. Despite 
that, the binary inclination range we found was fairly narrow (12°) and we were able to 
provide a new constraint on the mass of the bla ck hole in XTE Jl 1 18+480. Interestingly, 
the binary inclination of 68° ± 2° determined by IGelino et al.l ( 120061 ) is at the low edge of 
our inclination range, suggesting that they may have acquired their light curve data when 
the NIR non-stellar contribution was small. 



In an extensive study of the X-ray binary A0620-00, ICantrell et al.l ( 120081 . 120101 ) showed 
that the system existed in three distinct optical states even in quiescence. The authors 
found that a correct determination of the inclination relies on identifying the state of the 
system. This was found to be more important than the particular waveband where the 
measurements are made since non-stellar sources were present at both optical and NIR 
wavelengths. XTE J1118+480 also exhib its changes in the s hap e of its quiescent ellipsoida l 
light curve (e.g., comparing the data of IGelino et al.l 120061 and iMikolajewska et al.ll2005l ). 
The asymme try we observed in t he H- band light curve of XTE J1118+480 is similar to 
that seen by IMikolajewska et al.l (120051 ) in their J-band data. Significant changes in the 
quiescent infrared light curves of the black hole X-ray binar y GR O J0422+32 were al s o seen 
in two independent observations by iReynolds et al.l (120071 ) and IGelino fc Harrison! (120031 ) 
even though the mean K-band magnitudes did not vary significantly. From these studies, is 
clear that X-ray binary systems continue to harbor extensive activity and multiple accretion 
and/or outflow states even when fully in quiescence. 



5. Conclusions 

We have obtained broadband near-infrared spectroscopy of XTE J1118+480 and con- 
temporaneous light curve data to accurately account for the veiling that affects determina- 
tions of the binary inclination and compact object mass. By comparing the shape of the 
spectral energy distribution of the combined NIR spectrum as well as individual absorption 
lines in our spectrum with those of field stars of known spectral type, we were able to broadly 
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account for the fraction of NIR light contributed by the donor star as / = 0.50 ± 0.32. We 
factored the non-stellar contribution into our H-band light curve fits and obtained a binary 
inclination of 68° < i < 79°. From these we obtained a robust determination of the black 
hole mass in XTE J1118+480 of 6.9 < (M BH /M Q ) < 8.2. Our results are consistent with 
the picture of continued activity in X-ray binary systems even in quiescence and underscore 
the importance of accounting for all emission sources in the binary system when modeling 
spectra and light curves, even at NIR wavelengths. 

We would like to thank Emma Hogan for her assistance with Gemini-IRAF and Bernadette 
Rogers for letting us use her code to make GNIRS data compatible for use with Spextool. 

Facilities: Gemini-N (GNIRS), APO (NICFPS) 



A. Determination of error in the donor fraction calculation 

Due to the low SNR of our spectrum, the donor fraction obtained from the spectrum 
of XTE J1118+480 is only meaningful if accompanied by reliable uncertainty estimates. In 
order to adequately account for the noise in the spectrum, we created "scrambled" versions 
of the observed spectrum and used fits to the noise spectra to determine the uncertainties on 
the donor star fraction. For each of the J-, H-, and K-band spectra, we obtained the Fourier 
transform of the normalized spectrum. The Fourier-transformed spectra give a complex 

N 

valued array , X(k) = x(n)e^~ 27ri ^ N ^ n ^ k , where N represents the total data points, k 

n=l 

represents the individual Fourier transformed frequency components, and x(n) represents 
the normalized flux values. The amplitude of the Fourier transform is given by A{k) = 
\/ Re[X(k)] 2 + Im[X(k)] 2 and the phase is given by 4>{k) = tan" 1 (— ^^(h)) ) • 

To estimate the noise in our spectral data, we "scrambled" the phase of the Fourier- 
transformed spectrum by adding a different random number between and 2ir to the original 
phase at each k while retaining the original power spectrum |Xfc| 2 . We used the rand function 
in MATLAB that generates uniformly distributed random numbers between and 1. We 
then reconstructed the spectrum by calculating the inverse Fourier transform of the "phase 
scrambled" Fourier-transformed spectrum. In Figure [TTj we depict an example of calculating 
the donor fraction when an MO. 5 V star was compared to the reconstructed 'phase scrambled' 
spectrum of XTE Jll 18+480 in the H-band. 

The advantage of scrambling the phases in Fourier space is that, unlike scrambling 
the original wavelength-binned data, it works for both white and non-white noise. Any fit 
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done to the "phase scrambled" spectrum will be spurious since we are fitting signal to pure 
noise which allows us to place uncertainties on our fits to the original spectrum. We fit the 
spectra of field stars with K5 V, K7 V, MO/MO. 5 V and Ml V spectral types to each of 
the phase scrambled spectra and calculated the donor fraction for each wavelength region 
under investigation listed in Table [2j For each field star and each wavelength region under 
investigation, we repeated this process for 50 randomly generated phase scrambled spectra 
and examined the variance of the donor fractions. The square root of the variance is then 
a robust indication of the uncertainty in the donor star fraction. The errors listed with 
the donor fractions in Table [2] was obtained from this analysis. Finally, these errors were 
propagated through the remainder of the steps leading up to the average donor fractions in 
the J-, H- and K-bands. 
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Table 1. Observations of XTE Jll 18+480 





Date 


Individual t exp 
(sec) 


Total t obs 
(hr) 


Orbital Phase Coverage 


Spectroscopy 


04/02/11 


310 


2.1 


0.13 - 0.39, 0.72 - 0.98 




04/12/11 


310 


2.1 


0.06 - 0.32, 0.49 - 0.76 


Photometry 


04/02/11 


20 


4.0 


0.00 - 0.35, 0.40 - 0.99 




04/03/11 


20 


3.4 a 


0.18 - 0.99, 0.01 - 0.04 



a Data after the two hours of observation was discarded owing to bad weather. 
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Table 2. Donor star fraction fits 



Band 


Wavelength 


rangc(/xm) 


Dominant Feature 


/ (K5 V) 


/ (K7 V) 


/ (MO. 5 V) 


/ (Ml V) 


K 


2.203- 


2.214 


Na I 


0.80±.15 


0.70±.15 


0.65±.18 


0.62±.18 


H 


1.48 - 


1.51 


Mgl 


0.37±.07 


0.37±.07 


0.48±.08 


0.58±.09 




1.568 - 


- 1.60 


Mg I, Si I 


0.40±.ll 


0.44±.09 


0.47±.09 


0.63±.10 




1.67- 


1.68 


Al I 


0.39 ±.12 


0.39±.ll 


0.30±.09 


0.40±.10 




1.70 - 


1.72 


Mgl 


0.55±.09 


0.60±.09 


0.65±.ll 


0.70±.09 


J 


1.176- 


1.202 


Mg I, Fe I, Si I 


0.63±.10 


0.65±.13 


0.68±.13 


0.69±.15 




1.310- 


1.317 


Al I 


0.65±.15 


0.64±.16 


0.66±.17 


0.72±.18 



Table 3. Adopted Physical Parameters for XTE J1118+480 



Parameter Value Reference 



P orb 0.16995±0.00012 d Gonzalez Hernandez et al. (2008) 

K 2 708.8±1.4 km s" 1 Ibid, 

q 0.024±0.009 Calvelo et al. (2009) 

i 68° - 79° This work. 

M B h 6.9 - 8.2 M This work. 



Donor spectral type K7 V - Ml V This work. 
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Fig. 1. — The time-averaged spectrum of XTE J1118+480 obtained after correcting for 
atmospheric absorption and shifting the individual exposures to the rest frame of the donor 
star. The dereddened spectrum is shown, the SpeX short cross-dispersed observing mode 
does not cover the wavelength region from 1.86 - 1.93 /im, which is indicated by a straight 
line in the figure. 
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Fig. 2. — The J-band spectrum of XTE Jll 18+480 showing emission features from the 
accretion disk and absorption features from the donor star. The typical error bar is shown 
in red. 
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Fig. 3.— The H-band spectrum of XTE J1118+480. 
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Fig. 4.— The K-band spectrum of XTE J1118+480. 
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Fig. 5. — The time-averaged spectrum of XTE J1118+480 compared to the spectrum of a 
K5 V star (in red) and a Ml V star (in green). The K5 V and Ml V spectra have been 
scaled to match the XTE J1118+480 flux near 2.23 /zm. 
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Fig. 6. — The normalized H-band spectrum of XTE Jl 118+480 is shown in the upper panel 
with the normalized spectrum of a K7 V star overplotted (in red). The lower panel shows the 
same spectra after the K7 V star has been normalized to the best fit fractional contribution 
of 0.50. 
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Fig. 8. — The light curve of XTE J1118+480 after binning the observed data in phase bins 
of 0.03. The error bars are derived from scatter about the mean in each bin. For phase bins 
that contained fewer than 3 points, we used the largest error bar in the rest of the bins as 
the uncertainty. 
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Fig. 9. — The best fit light curve obtained by modeling a donor star (fit range: <j)= 0.0-0.5) 
with a constant extra flux component. The inclination is 68°. 
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Fig. 10. — Light curve models consisting of a donor star and an accretion disk along with 
a bright spot are fit to the full light curve of XTE J1118+480. The best fit parameters for 
this model are % = 78°, T disk = 3000^, T spot = 12000^, <j) spot = 85°, A<f) spot = 5° 
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Fig. 11. — An example of fitting an MO. 5 V stellar spectrum to the inverse "phase scrambled" 
Fourier-transformed spectrum of XTE Jl 118+480. The donor fraction variance obtained 
from fitting many such randomly phase-scrambled spectra were used to place uncertainties 
on the donor star fraction. 



